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Abstract 

Molecular dynamics calculations and optical spectroscopy measurements of 
weakly active infrared modes are reported. The results are qualitatively un- 
derstood in terms of the "motional diminishing" of IR lines, a process analo- 
gous to the motional narrowing of a nuclear magnetic resonance (NMR) signal. 
In molecular solids or liquids where the appropriate intramolecular resonances 
are observable, motional diminishing can be used to study the fluctuations of 
the intermolecular interactions having time scales of lpsec to lOOpsec. 
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Group theory of molecular symmetries is extensively used in infrared and Raman spec- 
troscopies to assign vibrational modes to point group representations. One of the simplest 
and most important arguments applies to molecules with a center of inversion symmetry: 
the vibrations can be divided into symmetric (gerade, denoted g) and asymmetric (unger- 
ade, denoted u) modes; only a u mode may have IR-activity and only a g mode is allowed 
in the Raman spectrum. Except for the simplest molecules, further analysis is needed to 
determine if the modes are indeed allowed; the u modes with other than vector symmetry 
are IR-forbidden, and the g modes with other than tensor symmetry are Raman forbidden. 

The symmetry group of an isolated molecule is strictly relevant only for low density gas 
phase samples. Intermolecular interactions break the pure molecular symmetry. In the gas 
phase, "collision induced scattering" (CIS) has been observed for some time, leading to an 
enhanced Rayleigh tail in noble gases |l|J^], a broadening of the Raman line in SF 6 || and 
the appearance of "forbidden" IR modes in the Raman spectrum of CO2 ||. In liquids the 
intermolecular forces are even stronger and a large body of literature discusses the various 
manifestations of these interactions |§-|[|. In the solid state the static perturbation leads to 
the well known crystal field splittings of degenerate modes, in addition to the appearance of 
previously "forbidden" resonances in the spectrum |TTj-|r3||. 

For disordered systems the treatment often starts in terms of two experimental time 
scales: the energy relaxation time 7\ and the dephasing time T 2 [[f|,|l5|. These times are 
analogous to T\ and T 2 in nuclear magnetic resonance (NMR) ||16|| . T\ tells us how fast 
the energy of a given vibration is dissipated to other modes; T 2 is the time required for the 
oscillators in the ensemble to develop a 2tt phase difference. Since in many systems Ti ^> T 2 , 
the width of a resonance line is, in good approximation, Au = 1/T 2 . 

The theory of the gas and the solid phases are relatively easy to handle, whereas the 
liquid presents a greater challenge. For CIS in the gas phase the dephasing is viewed as a 
result of independent collision processes, separated by the average time between collisions r c . 
Each collision adds a small random phase to the oscillations; consequently, Auo ~ 1/V C . The 
magnitude of phase shift depends on the "violence" of the collision relative to the strength 
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of the oscillator. A detailed calculation results in Auj = {k^T) / \uj 2 L 2 mT c ) , where k^T is the 
thermal energy, L is the range of interaction and m is a mass factor ||17|| . In the crystalline 
solid state the problem is formulated in terms of the exchange coupling energy V between 
the oscillators: the line splittings will be of the order of 5uj (xV. In a disordered solid the 
a linewidth Acu oc< V 2 > x l 2 is expected, where <> represents an appropriate average over 
all configurations. 

Approaching from these two extremes, the liquid can be either viewed as a dense gas, or as 
a disordered, fluctuating solid. A great amount effort has been concentrated on adapting the 
gas phase results to liquids Ul7j , including calculations taking into account the microscopic 
details of the interaction potential PJUfl. In this Communication we will take the other 
route. We will investigate what happens to the linewidth of the allowed resonance lines, and 
to the width and strength of the "forbidden" lines, if fluctuations are introduced to a quasi- 
static, random system. We perform computer simulations on a classical vibrational model 
and we use CS2 as a model system in the experiments. The agreement we find between the 
measurement and the calculations should not be viewed as an exhaustive and quantitative 
explanation for the behavior of the spectral lines, however the results do point to a novel 
contributing mechanism, which we term "motional diminishing" . 

Qualitatively, as the correlation time of the fluctuations (r) decreases, one expects that 
the effective symmetry breaking field will disappear. This is evident in the limit of extremely 
fast fluctuations, when r < T (where T = 2tt/uj, the period of oscillation of the molecular 
vibration in question) the external fields are expected to completely average out, producing 
effectively no symmetry breaking. Consequently, i. the linewidth of the allowed resonances 
narrows, Auj oc r []T5| | , and ii. the weakly allowed lines disappear. The question is, what 



is the critical timescale of the fluctuations for this crossover to occur? While the first 
effect ("motional narrowing") has been discussed before |T7,|I5|, the crossover from weak 



optical activity to no activity as the correlation time of the fluctuations decreases is a novel 
phenomena; we call it motional diminishing. 

In addition to shedding a new light on the behavior of molecular liquids, our consid- 



erations are directly applicable to "plastic" crystals [pl|. In fact the main motivation of 



this work originated from the desire to understand the behavior of a new and particularly 
interesting plastic crystal, Ceo [2D|-2^| . In a plastic crystal the molecules perform rotational 



diffusion, while preserving their center of mass positions at the crystallographic coordinates. 
If viewed from the reference frame of a given molecule, the magnitude of the crystal field 
does not change much but its direction seems to fluctuate with a correlation time r, the ro- 
tational diffusion time constant. Note that r and r c from CIS have opposite effects: Au oc r 
and Au; oc l/r c 

We performed computer simulations by solving the classical equations of motion for a 
simple, one dimensional, harmonic model system, in the presence of a harmonic driving field 
E of angular frequency u, and fluctuating external perturbations. The time averaged power 
dissipation was calculated by integrating the field x charge x velocity at various values of the 
drive frequency. The "ball and spring" system is similar to a CS2 or CO2 molecule (except 
that it does not model the bending mode): it consists of two balls of equal mass m and charge 
q attached to a third ball of mass M and charge Q by equivalent springs of spring constant 
k. This "molecule" has a center of inversion symmetry with one Raman-active g mode and 
one IR-active u mode. The external symmetry breaking perturbation is represented by two 
springs, k' and k", connected to the two outer masses. The Hamiltonian of this system is, 

tf = £ + i£ + lfr + - * 3 ) 2 + \Kx 2 - x 3 ) 2 

+ E[q(xi + x 2 ) + Qx 3 ] + \k'x\ + \k"x\ . (1) 

Randomness of the external perturbation is achieved by making the strengths of these outer 
springs fluctuate randomly. The fluctuations are characterized by a correlation time r; the 
probability of switching to a new random value of k' (or k") in time interval dt is dp =dt/r. 
All other effects of the environment, the coupling to other modes, etc., are represented by 
a damping force proportional to the velocity of the masses, with damping constant 7 for 
all three balls. This damping provides a finite width for each resonance mode even if no 
outside perturbation is provided. Typical numerical values were k — 1, m = M = 1, q — —1, 



Q = +2 and 7 = 0.005. 

If one side of the system is perturbed by an time-independent potential (represented by 
k'), then two resonance lines are observed: a weak one close to the "forbidden" g mode 
frequency to — ujq — ^Jk/m and a strong one around the u mode, to = \/3lj . The line 
shifts from the "ideal" positions are linear in k = k'/k; the line intensity of the "forbidden" 
resonance is proportional to k 2 . 

First we consider random fluctuations where a perturbing spring of fixed magnitude 
is randomly switched between the left and right sides, the unperturbed side having zero 
external coupling. In this system (referred to as "case I" later), the resonance mode frequency 
is the same for both the "left" and "right" perturbations, and therefore there is no additional 
line broadening due to the external forces. We find that slowly fluctuating perturbations 
produce a weak IR mode close to u = u , but this "forbidden" resonance vanishes as the 
correlation time r becomes less than the lifetime of the resonance r < Ti = I/7. Note that 
the line narrows even if r is much longer than the period of the oscillation. 

A more realistic model should involve random fluctuations in the magnitude of the ex- 
ternal perturbations (this will be case II). Figures [I] and |2] summarize the results obtained 
when the k! and k" spring constants were assumed to vary independently and randomly 
with a uniform probability distribution in the range of < k < K max = 0.4. For quasi- 
static perturbations (long correlation times such as r = 5000) the u mode and the g mode 
resonances are both broadened beyond the lifetime effects seen in case I (the quasi-static 
linewidth Acj(t — > 00) = 5 is larger than 7), and the center frequencies are shifted from 
the unperturbed frequencies to somewhat higher values. At the other time extreme, fast 
fluctuations (like r = 5) lead to a vanishing infrared activity of the g mode (Fig. |]). The u 
mode narrows to the "lifetime" width 7, and its integrated intensity remains approximately 
unchanged (see the upper inset of Figure Q). 

The parameters in this work are deliberately selected so that the three time scales in the 
problem (T = 2tt/uj, T* = 1/5, and T\ = I/7) are well separated. For the g mode, in the 
dimensionless units of the calculation, T = 1.102, T* = 18 and T\ = 200; for the u mode the 
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values are 1.758, 40 and 200, respectively. The insets in Figures |] and|l| clearly indicate that 
the crossover from the quasi-static behavior at long r (broadened IR line and non-zero IR 
intensity of the Raman line) to the fast fluctuations regime (narrowed IR line and vanishing 
IR-activity of the Raman line) occurs at around r = T* for each mode. The intensity of 
the "forbidden" Raman line as a function of r is well fit by the function I = I /(l + —), as 
demonstrated by the solid line in the inset to Figure [l]. In contrast to case I, the lifetime of 
the resonance, 7\, does not play a particularly important role in the physically more realistic 
model of case II. 

According to these results, the motional narrowing of the IR line is indeed very similar 
|7|JT5[| to the motional narrowing of NMR resonances [|l^,^]. The motional diminishing of 



the IR-activity of the g mode is, however, an entirely new phenomenon. We found that 
the inverse linewidth of the quasi-static perturbed resonance determines the crossover time. 
In case I this linewidth is due to lifetime broadening, in case II it is determined by the 
magnitude of the perturbation itself. Although our discussion was conducted in terms of 
a classical system, experience with quantum oscillators suggests that the main features of 
this crossover behavior would survive a full quantum mechanical treatment. All of these 
phenomena can be treated in a unified way in terms of the fluctuation-dissipation theorem 

We choose CS2 as a model system to test the predictions of our theory. In this material 
the "forbidden" Raman lines have been studied previously in detail as a function of tempera- 
ture P JT8|j24| and hydrostatic pressure [p5| -|28|| . The rotational diffusion of the molecules has 



been studied independently by NMR, and the correlation time is well known as a function 
of temperature p9|j30|| . 

We measured the IR transmission by using a 0.1mm thick liquid cell, enclosed by 15mil 
thick, singly polished, high-resistivity Si windows. The cell was filled with Fischer Scientific 
spectrophotometric grade CS2 liquid, sealed with teflon washers and mounted on the cold 
finger of a Helitran liquid N 2 flow refrigerator. The spectra were recorded in the 125K-300K 
temperature range by a Bomem MB-155 FT-IR spectrometer, using the spectra of the empty 



cell as backgrounds. As shown previously, the transmission method is well suited for the 
investigation of weakly-active modes ||31|| . 

The spectra around the v\ (Raman-active stretching) mode are shown in Figure |3|. At 
frequencies extending up to 3000cm -1 we clearly observe several resonances corresponding 
to combination modes (z/ 3 — u 1 , v x + u 2 , Ui + 3u 2 , Ui + z/ 3 , 2u 2 + z/ 3 , 2ui + z/ 3 , and Ui + 
2v 2 + ^3)- With the sample thickness used in this measurement the IR-active v 2 and v% 



modes appear totally saturated. Symmetry considerations |3^j33|] on the unperturbed, non- 
harmonic system predict that the combinations composed of an odd number of IR modes are 
IR-active. The intensity of these modes depends on the Bose factors rii = n(Tioji/k-QT) of the 
contributing fundamentals (i is the mode index of the fundamental). With the sole exception 
of the z/ 3 — Vi difference mode seen in Fig. |3], the intensities of the combination modes are 
expected to be, and found to be, approximately independent of the temperature. In Figure 
[| the integrated intensity of this difference mode is shown as a function of temperature 
(circles), together with the expected \n\ — n 3 \ temperature dependence (dashed line). 

The V\ and the 2v 2 resonances are IR-forbidden by the group theory (they are in fact 
the two strongest Raman modes). These lines appear in the IR spectrum either due to 
a symmetry-breaking from the interaction between the CS2 molecules or between CS2 
molecules and impurities in the liquid. At low temperatures both "forbidden" lines are 



clearly observed. The v\ resonance has a Fano lineshape [34J]. As the temperature is in- 
creased the lines become broader and their intensities decrease, in sharp contrast to the 
behavior of the "allowed" modes (independent of temperature) and the "difference" mode 
(increasing with temperature). However, inserting the known NMR correlation times into 
the motional diminishing results described above gives a complete account of the tempera- 
ture dependence (Fig. |], solid line). We emphasize that, other than matching the measured 
and calculated intensities at one temperature, there are no fitting parameters in the calcu- 
lated lines of Fig. |. 

In a similar manner we can interpret the results of a previous pressure dependence study 
of the Raman spectrum of CS2 by Ikawa and Whalley |25|-p7|. Increasing pressure in CS2 is 



analogous to decreasing temperature: the fluctuations slow down (r increases). The Raman 
observation of the IR-active v% mode being sharp at high pressures, and broadening and 
reducing in intensity at lower pressures (motional diminishing) can therefore be explained. 
Similarly the Raman-active polarized spectrum of the v\ mode is motionally narrowed, while 
the depolarized V\ is motionally diminished at lower pressures. 

To summarize, using the results of a computer simulated model, we qualitatively ex- 
plained the observed diminishing of IR-inactive vibrational modes in terms of a fluxuating 
molecular environment. While motional diminishing is not the sole mechanism which re- 
duces the intensity of weak spectroscopic features ||, we have shown that in appropriate 
situations its effect can be a significant contributor. The results discussed in this work can 
be utilized to measure fluctuating time scales in solids, liquids and glassy systems. In a 
more complex system with IR- and Raman-silent modes, all weakly active fundamentals 
may exhibit similar motional diminishing given proper symmetry breaking fluctuations. By 
investigating the temperature dependence of the intensity of these weak modes, molecu- 
lar fluctuations falling into the timescale range determined by the experimentally resolved 
linewidths of spectrometers, typically 300GHz-30THz, can be studied. 

This work was supported by NSF grant DMR9202528. 
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FIGURES 

FIG. 1. Time averaged power absorbed as a function of frequency, u and random fluctuation 
time, r, in the frequency range near the g mode. Inset shows the integrated intensity, its fit to the 
given functional form, and width of the peak as a function of r. 

FIG. 2. Time averaged power absorbed as a function of frequency, to and random fluctuation 
time, r, in the frequency range near the u mode. Upper inset shows the integrated intensity and 
width of the peak as a function of r. Lower inset displays how the character of the resonance 
changes from Lorentzian at low r to a mix of Lorentzian and Gaussian a high r. 

FIG. 3. A portion of the IR transmission spectrum of CS2 as a function of temperature. Labels 
denote the assignment of the lines, and the origin of the temperature dependencies. 

FIG. 4. Temperature dependencies of the integrated intensity of the v\ "forbidden" mode 
(squares) and the 1/3 — v\ difference mode (circles). Results from motional diminishing theory is 
shown as a solid line. The dashed line represents the Bose factor temperature dependence. 
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